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I. INTRODUCTION 

Freeze out (FO) is a term referring to the stage of ex¬ 
panding or exploding matter when its constituents (par¬ 
ticles) loose contact, collisions cease and the local dynam¬ 
ical equilibrium no longer can be maintained. When the 
local equilibrium is significantly perturbed, the micro¬ 
scopic length and time scales become comparable to the 
characteristic macroscopic ones and the hydrodynamical 
approach used to describe the evolution of matter breaks 
down. In the absence of collisions the momentum dis¬ 
tribution of the particles ’’freezes out”, hence the name 
kinetic freeze out. 

The final break-up corresponds to a ’’phase-transition” 
from an interacting fluid to a non-interacting gas of par¬ 
ticles, where the interactions between the constituents 
ceases suddenly when reaching the ” critical” FO temper¬ 
ature of the order of the pion mass Tpo ~ 140 MeV, 
as first assumed by Landau 0- Consequently, FO is a 
discontinuity in space-time represented by a space-time 
boundary or FO hypersurface, taken at the critical tem¬ 
perature [i.e., FO isotherm]. Across such FO hypersur¬ 
face the properties of matter change suddenly. We denote 
the two sides of the FO hypersurface as Pre FO and Post 
FO sides. Originally, the Post FO distribution function 
was assumed to be an equilibrated Jiittner distribution 
function boosted with the local flow velocity on the ac¬ 
tual side of the FO hypersurface. This approximation 
and method corresponds to the so-called ’’sudden” FO 
model by Cooper and Frye 0. 

The Cooper-Frye type of FO process is the zero thick¬ 
ness limit of a more realistic, so-called ’’gradual” FO pro¬ 
cess, where the FO description applies over a finite space- 
time domain, i.e. FO layer. Inside the finite FO layer the 
properties of the matter change gradually trough inter¬ 
actions, while the frozen out particles are formed and 
emitted at different ’’temperatures” which correspond to 
the actual temperature of the interacting matter, gradu¬ 


ally during the whole evolution of the matter. 

The basic philosophy of this paper is similar to the 
recent work [l[, which introduces and analyzes in detail 
the gradual FO description for space-like FO situations. 
Through the paper we are going to use the notation from 
Ref. Q, recall the governing equations and the major re¬ 
sults for comparison. We have made this paper sole and 
complete and widely understandable without the need to 
consult our previous work where the original ideas were 
first introduced. 

Time-like discontinuities represent the overall sudden 
change in a finite volume where the events happen si¬ 
multaneously at causally disconnected points of the hy¬ 
persurface with a time-like normal vector. For example 
the assumption of instantaneous or isochronous FO [i.e., 
happening at a constant time in the center of mass sys¬ 
tem], belongs to this category. 

The aim of this paper is to present an analyze a simple 
gradual kinetic FO process with time-like normal vectors. 
In the first part of this study we introduce and general¬ 
ize the gradual kinetic FO treatment for a finite time-like 
FO layer in a fully covariant footing, while in the second 
part we analyze its outcome. 


II. FREEZE-OUT FROM A FINITE TIME-LIKE 
LAYER 

The basis of the gradual FO method is to separate the 
’’full” / = f{x,p) distribution function into still inter¬ 
acting and already frozen out parts, / = p -\- , and 

describe the evolution of both components in a self con¬ 
sistent way siiai- This can be achieved by intro¬ 
ducing the so-called escape rate, Vescis,p), used to drain 
particles, which no longer collide from the interacting 
component, f\ and to gradually build up the free com¬ 
ponent, f^. For the better understanding of the model 
we use Fig. [U and assume that the FO of particles starts 
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from the inside boundary of the FO layer, (thick line). 
Within the FO layer of finite thickness, L, the density of 
interacting particles decreases and disappears once we 
reach the outside boundary, S 2 (thin line), of the FO 
layer. 

In general, the kinetic description of freeze out leads to 
a complex multidimensional problem. To clarify the ba¬ 
sic properties of the FO process through a finite layer we 
may essentially reduce the number of variables, assum¬ 
ing that the dominant change in the distribution function 
happens in the direction of the FO normal vector, while it 
is negligible along the directions perpendicular to it, (e.g. 
in a spherically symmetric system the change happens in 
radial direction, and it is negligible in the perpendicular 
directions). Thus, the FO process can be effectively de¬ 
scribed as a one-dimensional process and the space-time 
domain where such a process takes place can be viewed 
as a FO layer, where the FO normal vector is tied to the 
direction of the density decrease arising from velocity di¬ 
vergence at a curved FO surface. 

If we have a space-like normal vector, da^ = (0,1, 0,0), 
the resulting equations can be transformed into a frame 
where the process is stationary, while in the case of a 
time-like normal vector, = (1,0, 0,0), the equations 
can be transformed into a frame where the process is uni¬ 
form and time-dependent. In this paper we only discuss 
FO processes inside a finite time-like layer. 

Here we recall the governing equations [i.e., Eqs. (13- 
14)] from Ref. [H, which can be used in both time-like 
and space-like FO cases. The equations depend on the 
projection, s = in the direction of the FO normal 

vector, da^, where the four vector denotes the particle 
coordinate, having its origin at the inner surface of the 
FO layer. Thus, 

dsf{s,p) = -Vesc{s,p) f{s,p) + ^ ^ 

Mh 

dsfHs,?) = +'Pescis,p) f{s,p) , ( 1 ) 


where using the relaxation time approximation we ensure 
that the interacting component approaches the equili¬ 
brated Jiittner distribution, /*g, with Xth relaxation time 
(or relaxation length in the space-like case). 

The escape rate, Vescis,p), describes the escape of par¬ 
ticles from the interacting component into the free com¬ 
ponent and it is defined as: 


esc (-^j P') 


1 

A(s) 


L p^da^ 
L - s pt^Uf^ 


<d{p^da^) , 


( 2 ) 


where the parameter, L, is the ’’proper” thickness of the 
FO layer and it is an invariant scalar. The proper thick¬ 
ness is analogous to the proper time, that is the time 
measured in the rest frame of the particle. In our case 
the local rest frame is the rest frame attached to the FO 
front (RFF), (see III A|) . and thus the proper ’’thickness” 
of the FO layer is the invariant proper time interval be¬ 
tween the start of the process and its end, (see Fig. [1] 
at point A). Furthermore, p^ is the four-momentum of 



FIG. 1: (Color online) The figure shows a finite FO layer 
with varying thickness. The normal to a surface element is 
dafj,, thus between A and B the surface is time-like, [i.e., 
da'^dafj, = -1-1], from B down to C it is space-like, [i.e., 
da'^dafj, = — 1]. The change from time-like to space-like sur¬ 
face happens where the normal of the surface is light-like, 
but not necessarily has its origin at the center of the system. 
The momentum of particles is and the four vector de¬ 
notes the particle coordinate, having its origin at the inner 
surface of the FO layer. On the time-like FO region all parti¬ 
cles emerging from a point on the Pre FO side will propagate 
to the Post FO side, while the particles originating from the 
space-like part of the FO surface are divided between Pre FO 
and Post FO parts. Only those particles cross the surface 
which have their momentum enclosed by the light cone and 
the Post FO surface, [i.e., if p^da^, > 0]. 


particles, da^ is the normal in the FO direction, while 
is the flow velocity normalized to unity. The initial 
characteristic time is denoted by tq, and the ©(p^du^) 
function, was first introduced by Bugaev @ to ensure 
that all particles leave to the outside. In case of time-like 
FO this condition is always satisfied and does not lead to 
any additional constraint. 

A qualitative expression of the escape rate for both 
time-like and space-like FO situations is based on the 
following simple assumptions. Particles with higher mo¬ 
mentum in the FO direction will freeze out first. The 
particles closer to the outside boundary of the FO layer 
have a greater chance to freeze out since the probability 
to find another particle to collide with is smaller as the 
system became sparser, (given by the L/{L — s) factor). 
In our special case the FO direction is parallel to the time- 
axis, thus all particles will freeze out irrespective of their 
momenta. However, particles emitted at later times will 
freeze out ’’faster” since they have less chance to collide 
with other particles in the diluted system. Please note 
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that, although tq is assumed to be a constant for simplic¬ 
ity, the characteristic FO length is increasing with time 
or distance such as, tq{L — s)/L. The detailed treatment 
and analysis of the escape rate can be found in [iiMini. 

Now, if we describe the time evolution of the particle 
FO, then dcr^ = (1,0,0, 0), x^da^ = t and p^da^ = 
thus eq. Ill) leads to 


dtp 

dtff 


1 

To 

1 

H- 

To 


L-t 

L 

L-t 
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|0 




P^Ufj, 


r + ^ 


fi, - f 


'^th 


(3) 


where the interacting component approaches the equili¬ 
brated Jiittner distribution, /eg, with Tth relaxation time. 
This is a common simplification and the practical reason 
to use it is to calculate the quantities depending on the 
equilibrium distribution function. Although the present 
solution mathematically is achieved taking an infinitely 
short relaxation time, in reality, one can show, see Ref. 
[l[, that the complete thermalization of the interacting 
component can be achieved with good accuracy if tq is 
smaller than r by a factor of 2 or more. For a thorough 
analysis of this approach, see Refs. and the Ap¬ 

pendix. Here we mention that in our calculations we use 
only one type of particles, namely massless pions, there¬ 
fore the chemical composition of our system remains un¬ 
changed during the kinetic FO. Furthermore, for simplic¬ 
ity we assume simultaneous chemical and thermal equili¬ 
bration, thus the chemical potential of the massless pion 
gas is /X = 0 during the FO process. 

Earlier in Ref. [T^, the gradual time-like FO descrip¬ 
tion was modeled with equations having a similar form 
to eqs. ©, but it was only treating the simplest case 
when = dcr^, while the FO was lasting infinitely long. 
The model was based on the idea of the boost invariant 
Bjorken hydrodynamical model [l^ . where the evolution 
of matter is a function of the proper time, r, only, while 
the flow of matter is parallel to the normal vector of the 
proper time hyperbolas in every point. The covariant 
equations, eqs. m, return those equations if we change 
the time variable, dt = dr, and neglect the L/{L — t) 
factor. However, the new equations allow ^ da^ and 
also a possibility to finish the FO process within a given 
duration. 



FIG. 2: (Color online) A simple FO hypersurface in RFG with 
coordinates [t,x], where m'' = (1, 0, 0, 0)i{fg- The normal 
vectors of the FO front, da^, are time-like at points. A, B, 
C, while the normal vectors are space-like at points, D, E, 
F. At point B in RFF with coordinates [t’,x’], where da^ = 
(1,0, 0, 0)rff and = 7 (t(1, —Va,0, 0 )rff- Note that RFF 
moves together with the FO front, while on the figure the 
origin of RFF is shifted to match the origin of RFG. 


normal vector to the FO hypersurface for the time-like 
part is always da^ = (1,0, 0, 0 )_ rff , while on the space¬ 
like part is always = (0, 1 , 0 , 0 ) fff - 

If we are in the RFG then the four-flow is always 
= (1,0, 0, 0 ) ffg - If we take different characteristic 
points, for example points. A, B and C, on the FO hy¬ 
persurface then the normal vector is different at differ¬ 
ent points of the hypersurface in RFG, see Fig. To 
calculate the parameters of the normal vector, dtr^, for 
different cases in the RFF we make use of the Lorentz 
transformation. 

The normal vector of the time-like part of the FO hy¬ 
persurface may be defined as the local t'-axis, while the 
normal vector of the space-like part may be defined as 
the local x'-axis. This defines the axes of RFF. 


A. Reference frames 


B. Conservation laws 


Before proceeding further, first we define the reference 
frames in which our calculations will be handled. 

On the Pre FO side the matter is parameterized by an 
equilibrium distribution function, as required by the hy¬ 
drodynamical description of evolution. The frame where 
the matter is at rest is the local Rest Frame of the Gas 
(RFG), where = (1, 0,0, 0)ffg- On the Post FO side 
we can use the frame which is attached to the FO front, 
that is the Rest Frame of the Front (RFF). In RFF the 


The change of conserved quantities caused by the par¬ 
ticle transfer from the interacting matter into the free 
matter can be obtained in terms of distribution function 
of the interacting matter calculated from eqs. m as: 


dNtit) = dt[ ^p^dtP 


(4) 


dt L f d^p p^dan 


pPUp 


Penit.P). 


To L-t 
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while the change in the energy-momentum as: 


dT^it) = dtj 


dt L 
To L-t 


d^p 


p^p’' 


P^dap 

pPUp 


feq{t,P) 


(5) 


The equilibrium distribution function for massless bary- 
onfree particles is: 


fLit,p) = 


{2TrhY 


exp 


7 (p° — jup cos Op) 


T 


( 6 ) 


where the four-momentum of particles is p^ = (p°,p), 
p = IpI, = pcosdp, the flow velocity of the interacting 
matter in RFF is = 'y{l,v,0,0)iiFp, 7 = 1 /yr^, 
■*^=^ 11 ^ = sign(ti), and g is the degeneracy of particles. 
The results of the calculations in the RFF can be found 
in Appendix B. 

The change in energy density after a step dt is 


dei{t) = Upp{t) dT['''{t) , (7) 

from which by using a simple EoS, e = assT'^, for a 

2 

baryonfree massless gas with ctsb = ^, we can calculate 
the change in the temperature and Landau’s flow velocity 
similarly to [3,1,0. Thus, 


dlnT 

dv 


7 


4asBT4 

3 


dT/^ - 2vdTr + v^dTi 


- vdT^'^ -f (1 -h v^)dT]'^ - vdTi 


( 8 ) 


while in the massless limit the above equations lead to: 


dlnT 


dv 


dt / L \ 3n7 

To \L- t) AassT^ ’ 
dt / ^ \ 

To \L-t) A'-fassT^ ' 


(9) 


III. RESULTS AND DISCUSSIONS 

In this section we will present the results for the Post 
FO distribution and the relevant quantities calculated 
form this model. We will present our results for two 
different cases, for infinite FO (I) and finite FO (F). 

I) The system is characterized by an infinitely long 
FO duration, [i.e., L/{L — t) ^ 1 and {t = IOOtq), 
where most of the matter is frozen out]. The results 
are shown on Figs. [3l [5l [3 [H 

F) A finite FO process happening in a finite FO layer, 
where (L = IOtq). The results are shown on Figs. 

SUHElITol 



FIG. 3: (Color online) The temperature of the interacting 
component in RFF, calculated for an infinitely lasting FO. 
The initial temperature is To = 170 MeV, the parameter, vo, 
is the initial flow velocity. This corresponds to case I. 


A. The evolution of temperature of the interacting 
component 

The first set of figures, Fig. [3] and Fig. 13 shows the 
gradual decrease in temperature of the interacting com¬ 
ponent calculated in RFF. 

Comparing Fig. [3]with Fig. [3 we see the difference be¬ 
tween the finite and infinite FO. The FO in a finite layer 
is faster than in an infinite layer, ’’per se”. The temper¬ 
ature curves belonging to different initial flow velocities 
but with opposite sign [i.e., vg = —0.5 and vg = 0.5] are 
the same. This is so since for time-like FO we do not have 
any constraint on the momenta, such as the cut-off factor 
Q{pPd(jp), and the initial momentum distribution is sym¬ 
metric over the time axis. In the case of time-like FO the 
gradual cooling of the matter is faster and ’’smoother” 
compared to space-like FO, since the most energetic par¬ 
ticles freeze out unrestricted in direction thus the remain¬ 
ing interacting component cools down faster. 

The matter with higher initial flow velocity, vg, cools 
faster, but for small differences between the initial flow 
velocities, the resulting difference in the temperature is 
negligible. If the interacting gas has a higher flow veloc¬ 
ity in RFF, the escape rate is bigger for higher values 
of the flow velocity. This expresses the fact that if the 
matter flows we remove energy faster form the interact¬ 
ing component. 


B. The evolution of common flow velocity of the 
interacting component in RFF and RFG 

The second set of figures. Figs. [3 [3 shows the evo¬ 
lution of the flow velocity of the interacting component 
calculated for a baryonfree massless gas in RFF. 

Again, comparing Fig. [3 with Fig. [3 we can see the 
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FIG. 4: (Color online) The temperature of the interacting 
component in RFF, calculated for a finite {L = lOro) FO 
time, where the initial temperature is To = 170 MeV, and vo 
is the initial flow velocity. This corresponds to case F. 



FIG. 5: (Color online) The evolution of the flow velocity of the 
interacting component calculated for an infinitely lasting FO, 
corresponding to case I. The initial temperature is To = 170 
MeV, and vo is the initial flow velocity of the gas. 


difference between finite and infinite FO. In the case of 
finite FO the velocity decrease is much faster than in the 
case of infinite FO. 

Furthermore, we notice that the flow velocity of the 
interacting component tends to zero, while here we recall 
to compare, that in case of space-like FO it tends to —1. 
Again, this is due to the cut-off factor which retains par¬ 
ticles propagating with negative momenta in space-like 
directions. In RFF the quantities change discontinuously 
at the light cone and that is why we have different results 
for the final flow velocity comparing space-like and time¬ 
like cases. However, in RFG all quantities are continuous 
when crossing the light cone, see Ref. [10. 



FIG. 6: (Golor online) The evolution of the flow velocity of the 
interacting component calculated for a finite (L = lOro) FO, 
corresponding to case F. The initial temperature is To = 170 
MeV, and vo is the initial flow velocity of the gas. 


C. The transverse momentum and the contour 
plots of the Post FO distribution 

The third set of figures, Figs. [71 and [51 shows the evo¬ 
lution of the local transverse momentum distribution and 
the corresponding contour plots of the Post FO momen¬ 
tum distribution. 

We have presented a one-dimensional model here, but 
we assume that it is applicable for the direction trans¬ 
verse to the beam in heavy ion experiments. The plots 
presented should be related to the transverse momentum 
distribution of measured particles. 



FIG. 7: (Color online) The local transverse momentum (here 
Px) distribution for a baryonfree and massless gas at {py = 
0). The calculations were done for an infinite FO with thin 
lines and for a finite FO {L — lOro) with marker lines. The 
initial flow velocity and temperature are, — {1, 0 , 0, 0)rfg 
and To = 170 MeV. The transverse momentum spectrum at 
the end is obviously curved due to the FO process for low 
momenta. 

From Fig. |7|we see that the Post FO momentum dis- 
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FIG. 8: (Color online) The Post FO distribution, ffree{x,p), 
at point A of Fig. [2l for an infinitely long FO length. The fig¬ 
ures correspond to different time points, t = Itq, lOro, IOOtq 
respectively. Contour lines are given at values represented 
on the figure. The initial flow velocity and temperature are; 

= (1,0,0,0)flFG and To = 170 MeV. The maximum is 
increasing with t as indicated on Fig. 0 


tributions for the infinite and finite FO cases are quali¬ 
tatively identical. At the early stages of the FO process 
(for values of t ~ tq) the distribution of particles in the 
two cases match. This property persists until the end 
of the FO process, and can be seen on Fig. [71 where 
the local transverse momentum distribution calculated 
for an infinite and finite FO are identical. The maximum 
is increasing with t as indicated in Figs. |7]and|51 Thus, 
the final Post FO distributions do not differ if we switch 
from an infinitely long to a finite layer FO description for 
any initial flow velocity. This means that our finite layer 
FO description was done correctly. These important fea¬ 
tures of the model were already discussed in detail in 
Refs. [Il,[i|. 

The overall conclusion is that the resulting Post FO 
distributions are non-thermal distributions even for time¬ 
like FO processes. The distributions strongly deviate 
from thermal distributions in the low momentum region 
[i.e., px < 300 MeV]. If one decreases the duration of 
total FO time, the gradual FO process would still pro¬ 
duce similar particle spectra until the duration is not 
less than 2ro. Below that value the spectra becomes less 
curved and in the limit when the duration approaches 
zero the final momentum spectra corresponds to a con¬ 
stant temperature equilibrium distribution function. For 
L < 2ro the FO process does not have enough time to 
significantly change the shape of the final spectrum. 

Here one may go further and intuitively say that the 
FO process has a maximal lifespan, even though the pa¬ 
rameter, L, was not defined in this work in such way that 
it would allow us to exactly calculate its limits from first 
principles. Of course for such a statement to hold one 
would need a realistic full scale fluid dynamical simula¬ 
tion including chemistry, secondaries and the expansion 
of the system. However, our results concluded from this 
simple model could still hold valuable in the realistic FO 
modeling in complex fluid dynamical simulations. 
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FIG. 9: (Color online) The final Post FO distribution, 
ffreeix,p), at points A, B, C, of Fig. (2] calculated for a 
finite FO time, L = IOtq. The contour plots correspond to 
different initial flow velocities, = 7<t(1, —Va, 0, 0)rff, with 
Uct = 0,0.5c, 0.9c, respectively where the initial temperature 
is To = 170 MeV. Note that, the Post FO distributions at 
points B and C are not the boosted distributions of point A. 
The difference is that the Post FO distribution is much less 
elongated than the boosted Jfittner, because it is a superpo¬ 
sition of sources with decreasing speed in RFF as indicated 
in Fig. [6] 


D. The boosted Post FO distributions 

The fourth set of figures, Fig. |5| and Fig. (TUI shows 
the final Post FO distributions calculated for differ¬ 
ent flow velocities in RFF and the boosted post FO 
(Jiittner) distributions in RFG. The FO distributions 
corresponding to different initial flow velocities, = 
7ct(1, —Uct,0,0)/jff, generally lead to non-equilibrated 
and anisotropic Post FO distributions. However, boost¬ 
ing the distribution from point A to points B or C leads 
to a more elongated FO distribution in the direction of 
the boost than the calculated Post FO distributions at 
those points, (compare the contours with the same values 
given on Figs. l^and fTU)) . 

This is also an important outcome of our analysis lead¬ 
ing to the conclusion that assuming an isotropic equi¬ 
librated Jiittner distribution at time-like parts the FO 
hypersurface, other than at point A on Fig. where 
= da^, in general cannot hold [I^. More importantly, 
the common practice of boosting the Jiittner distribution 
or any post FO distribution function instead of calculat¬ 
ing it from the conservation laws, similarly as it was done 
here, leads to a noticeable difference in the final particle 
spectra. 


E. The non-equilibrated post FO distributions 
revised 

Here we present another important result of our study, 
following the approach from Ref. , where an infinitely 
long FO was studied with momentum independent es¬ 
cape rate. We can reproduce that earlier result by taking 
= (1, 0,0,0), and can calculate the temper¬ 
ature decrease using eqs. ®, in the case of a massless 








baryonfree matter for an infinitely lasting FO: 


V =0 


V = 0.5c 


V = 0.9c 


T(t) = T(to) exp - to)^ , (10) 

where k = 3/(47r^crsB). In the general case, the flow 
velocity is not zero, hence one has to solve the system 
of equations from eqs. ®. The distribution function of 
interacting particles (in the fast rethermalization limit) 
at any time t is: 

/■(t.P)=p^exp(-^e‘<-../™), (11) 

Now, we can solve the equation for the free component 
from eqs. @, therefore the distribution of free particles 
at time t is: 



P, [MeV] 


FIG. 10: (Color online) The post FO distribution at point 
A and boosted to the frames at points B and C as de¬ 
picted in Fig. [21 The different figures correspond to dif¬ 
ferent initial normal vectors, dcr'' = 7o-(l, Ua, 0, 0 )hfg, where 
Uct = 0,0.5c, 0.9c, and the initial temperature is To = 170 
MeV. 


■^0 Jtn 


( 12 ) 


Ei - 


k-^ 

- (2^3) L' 

which for t —>■ oo leads to: 
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P^Ufj,ito) 

Tito) 


(13) 


where Ei is the exponential integral function defined in 
Appendix B eq. (EOl)- Thus, we got a simple formula, 
similarly to the one in Ref. [l3| , which correctly param¬ 
eterizes the non-equilibrated post FO distribution func¬ 
tion when = (1,0,0, 0). It was actually shown in Ref. 
[ij that the post FO distribution is not sensitive to mo¬ 
mentum dependence of the escape rate, so we assume 
that this simple formula is valid for any initial flow ve¬ 
locity M^(to). 

Here we will use the above formula to plot the post FO 
distribution, for finite layers, with L > 2to, and extend 
this approximation to the general case when u^(to) 

On Fig. [TTl we have plotted the final FO distribu¬ 
tion functions calculated using eq. 031 ) with lines, and 
the finite FO (L = 3 to) calculation with marker lines. 
The results are matching, which is a remarkable result, 
thus we conclude that this simple approximation is ap¬ 
plicable for the description of gradual FO thorough finite 
time-like layers and correctly approximates its post FO 
distribution functions. 


IV. CONCLUSIONS 

In this work we have presented a simple kinetic freeze 
out model for a finite time-like layer. We have demon¬ 
strated that FO across time-like surfaces leads to non- 
equilibrated and anisotropic distributions. These distri¬ 
butions in general cannot be Lorentz transformed to a 
frame where the distribution is isotropic. The only ex¬ 
ception is when the normal to the FO hypersurface is 



FIG. 11: (Color online) The local transverse momentum (here 
Px) distribution for a baryonfree and massless gas at {py = 0). 
The calculations were done for an inhnite FO with thin lines 
using eq. m and for a hnite FO {L = 3ro) with marker 
lines. The initial temperature was To = 170 MeV. The differ¬ 
ent initial flow velocities are given on the flgnre legend. 


parallel to the local flow velocity. Our analysis shows 
that the usual practice of assuming a Jiittner distribu¬ 
tion as a Post FO distribution is in general not valid! 

We can also see that while the boosted Jiittner dis¬ 
tribution is elongated in the boost direction, i.e. in the 
direction of the Post FO distribution is close to a 
spherical and isotropic distribution at low momenta, and 
becomes elongated only at higher momenta, see Fig. |9l 
This special Post FO distribution leads to a curved 
- spectrum”. Here^we can also demonstrate (as in ear¬ 
lier works that non-equilibrium processes in 

kinetic FO lead to observable effects. 

We observe that the Jiittner distribution is not a good 
approximation for the Post FO distribution, just like in 
the case of a space-like FO. While in the case of space¬ 
like FO the Cancelling-Jiittner distribution introduced in 
Ref. [13 is satisfactory, in the case of time-like FO, we 
have found a simple formula to use. 

Now, one may ask the question whether we observe this 



















additional low pt effect in the experimental data. This 
effect has several possible explanations: products of low 
momentum resonance decays, the transverse expansion 
of the system, and possibly due to the long gradual FO 
with rethermalization. As already discussed, during such 
scenario the particles are freezing out at different grad¬ 
ually decreasing temperatures, thus correspondingly the 
final FO spectrum is a superposition of thermal distri¬ 
butions with different temperatures. Although the low 
Pt enhanced non-thermal spectrum of massless pions is a 
necessary outcome of long gradual FO with rethermaliza¬ 
tion, for the heavy particles (if these are in the mixture 
with pions) it is almost unobservable, see Ref. p^ . 

In our simplistic study we have found that FO in layer 
of finite thickness below L < 2ro will not show a sharp 
peak at low momentum in the transverse momentum 
spectrum. Thus, naively one can conclude from a simple 
fit that the FO in heavy-ion collisions happens in a nar¬ 
row or wide FO layer. However, such conclusion would be 
premature without including the expansion of the system 
and calculate two particle correlations. At the moment 
FO in a long finite layer, L > 2ro, cannot be excluded. 
If one assumes gradual FO with non-thermal post FO 
spectra, then one may also fit the data but with different 
flow velocity and slope parameter, where the curvature 
of the pion spectra at low pt will be partly due to FO 
and partly due to the expansion of the system. 


V. OUTLOOK 

We do not aim directly to apply the results presented 
here to experimental heavy ion collision data, instead our 
purpose was to study qualitatively the basic features of 
the freeze out process, and to demonstrate the applicabil¬ 
ity of this covariant formulation for FO in a finite layer. 

Here we note that our model may be applicable in CFD 
calculations, where one has both time-like and space-like 
parts of the full FO hypersurface, thus the gradual FO 
calculation must be done over the full FO hypersurface, 
with varying flow velocities and normal vectors. The 
method should be applied after reaching the Tpo crit¬ 
ical temperature and calculate the Post FO momentum 
distribution function starting form the inner FO hyper¬ 
surface. Such a calculation should be compared to the 
Cooper-Frye ansatz in the first place and then to exper¬ 
imental results, similarly to Refs. [l-flslfT^. 

A successful application of this model was already used 
to study the impact of nucleon mass shift on the freeze 
out process [l^ . This analysis will be carried forward to 
calculate the impact of mass shift on other particles, such 
as pions and kaons, which will help us study the effect 
of mass shift on two-particle correlations. An even more 
interesting study based on our analysis will estimate the 
effect of expansion on the time-like freeze out process us¬ 
ing the Bjorken model [2^. Therefore, we believe that 
our model may give a better description and understand¬ 
ing of the final observables which are calculated using the 


PreFO 


^th 

^f/j 


Post FO 


ao' 

1 


a ' 


X 


^0 


L = 4Ao 

FIG. 12: (Color online) A schematic view of the FO process 
for a linear density profile. Initially the mean free path is 
Am/p = “ 0 , the relaxation length is Xth = ‘iXmfp, the initial 
characteristic length is Ao, while the length of the FO layer is 
L = 4Ao. 

single (and two) particle distribution functions. 
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APPENDIX A 

Here we discuss the properties of the rethermalization 
term from eq. ([3]) and its consequences on finishing the 
FO process in a finite layer. 

From kinetic theory we know that if the following con¬ 
ditions: 

Un/p ^ '^th ^ 'h/ or ^mfp ^ ^th ^ Aq , (^4) 

between the average length between the collisions, the 
relaxation length, and the characteristic length are satis¬ 
fied, then we can use the Boltzmann Transport Equation 
(BTE) for the evolution of the single particle distribution 
function, f{x,p). 

For better understanding we first assume a linear de¬ 
crease of the interacting particle density during freeze 
out, such as, n(a:) = {L — x)/L, where the mean free 
path of interacting particles is, \mfp{x) « \/n{x). Using 
the relaxation time approximation for the FO process, for 
example at a; > 2Ao the density of the interacting par¬ 
ticles already decreased to n{x) < no/2, while the mean 
free path increased to Xmfp{x) > a = 2ao, see Fig. [121 
Consequently, by the end of the FO process, the thermal- 
ization length becomes longer than the initial character¬ 
istic length of the system. Although, tq is constant (scale 
parameter), the characteristic FO length is actually not 
constant during the evolution. Thus, using the rether¬ 
malization approximation the error we introduce within 
is of the order of Ttu/TQ. 
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In our model the change in the density is generally 
given as: 


dni{x) = dN^{x). 


(15) 


This leads to an exponentially fast decrease of particle 
density, therefore more than 95% of the interacting mat¬ 
ter is frozen out before Tth — tq, thus we can safely use 
the relaxation time approximation in our calculations. 


drnt) = - 


dT^^{t) dt L nT 


2 To L — t 8juj 
- Gt(m) + {m)\ 


m—O dt 1 
To 2 


dT^^it) +dT^^it) 


APPENDIX B 

The changes in the particle four current and energy 
momentum tensor are: 

dN°{t) = - GUm) + G^{m)\ 


dN^t) = 


To L-t Ajuj 

m—O dt L 

To L-t 

dN^jt) 

ju 

dt L n 
To L-t 4ju'y 

m—O dt L 
To L - t 


(3 -I- 


- 2b 


2Ki(a) + aKo{a) 


(3 + ?;2) 2 1 


5v 


(16) 


dT™(t) = - 


dt L nT j 


To L-t 4:ju-f ( 

m—O dt L 


- Gt (to) G^ (r 


-nT 


3(1-hu^) 7^ 


dTtit) = 


To L-t 
dT^°{t) dt L nT 


ju To L-t Aju'j 
- 262(3 -h u^)K 2 {a) - 2ab^Ki{a) 


m—O dt L 


To L-t 


nT 


3(1-1-1)2)73 7(3-I-1)2) 


dT^it) = 


dTf^{t) T 


JU 


lju 


dN^{t) - 


dN^jt) 

ju 


dt L nT f 252 „ 2 \ -r^ f \ 

- 7 — 

— 2juab^Ki{a)\ 


m—O dt L 


To L-t 


-nT 


3(1 -I-1)2)-)/3 


-h 


1 (l-l-3i)2)7 


71 )^ 


and 


dTr(t) = dTyy{t), 


(17) 


where a = ^, 6 = 07 and n = 'iTT^a?K 2 {a) g 
the particle density, while g is the degeneracy factor. Fur¬ 
thermore, the definition of the modified Bessel function 
of the second kind Kn{z) for n > — 1 , is 


Kr.{z) = T^,z- dxe-"(x 2 -z 2 )"-,. (18) 

( 2 n)! 


The analytically not integrable functions G~(m) and 
G+(to), where n > — 2, depend on other quantities such 
as u and T. However, in the case when, m —> 0, the de¬ 
pendence on other quantities persists and hence we only 
denote the mass dependence of the functions defined as: 


1 / _\ n 

Gt(m) = dppi^^/p"^ +m‘^j (19) 

The exponential integral function is defined as: 


fOO -X 

Ei(z) = / dx - 

Jz X 


( 20 ) 
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